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SUMMARY 

A rapidly convergent method is described for the calculation 
of the equilibrium composition of the carbon- hydrogen- oxygen-nitrogen 
system in the temperature range from 2000° to 5000° K. Charts are 
provided for the estimation and progressive adjustment of two inde- 
pendent variables on which the calculations are based. Additional 
charts Eire provided for the graphical calculation of the composi- 
tion. The use of the method and the charts is illustrated by an 
example . 


INTRODUCTION 

The determination of the equilibrium composition of a chemical 
reaction is necessary for the solution of certain problems in chem- 
istry and thermodynamics . When a reaction occurs at elevated tem- 
peratures, the products of reaction contain dissociated molecules 
which, when present in sufficient quantity, complicate the solution 
of the problem. The degree of dissociation is determined by the 
temperature of the products, their partial pressures, and an equi- 
librium constant for each dissociation reaction. 

For reactions involving carbon, hydrogen, oxygen, and nitrogen 
at temperatures between 2000° and 5000° K emd at pressures below 
100 atmospheres, a direct solution of the equilibrium equations is 
impractical if not impossible. Accordingly, some method of succes- 
sive approximation or trial and error is required. 

If only two elements Eire present in the reaction the problem 
is less difficult. For example, a direct method for computing the 
equilibrium composition of the hydrogen-oxygen reaction is described 
in reference 1. 
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A method of computing the equilibrium composition of a mixture 
of the elements, carbon, hydrogen, oxygen, and nitrogen, by means 
of a method of successive approximation is given in reference 2. 

The method vas applied to combustion in the cylinder of an internal- 
combustion engine. This method, although convenient for tempera- 
tures up to approximately 2000° to 2500° K, becomes cumbersome for 
temperatures at which a large fraction of the molecules are dis- 
sociated. 

A somewhat similar method for the same elements is described 
and the results are applied to the combustion in a gun in refer- 
ence 3. A compilation of thermal data required for equilibrium- 
composition calculations at temperatures up to 5000° K was given 
in this reference; a part of that data has been used in the pre- 
paration of the enclosed charts and has been reproduced as table I. 

A procedure is presented for computing the equilibrium composi- 
tion of a gas mixture corresponding to an assigned pressure and 
temperature that is especially convenient for the carbon- hydrogen- 
oxygen-nitrogen system in the temperature range from 2000° to 
5000° K and for pressures from approximately l/4 to 100 atmospheres. 
Charts that guide in the process of successive approximation and 
aid in the details of the calculations are presented. The use of 
the method developed at the NACA Cleveland laboratory is illustrated 
by an example. 

The products of reaction considered are gaseous water HgO, 
hydrogen Hg, hydroxyl radical OH, atomic hydrogen H, nitro- 
gen Ng, nitric oxide NO, atomic nitrogen N, carbon dioxide C0g, 
carbon monoxide CO, oxygen 0 2 , and atomic oxygen 0. For the 

temperatures, pressures, and chemical elements considered, the 
perfect gas laws were assumed and density correction for equilibrium 
constants was neglected. 


GENERAL SOLUTION 

The procedure for determining the equilibrium composition of 
a chemical reaction at a selected temperature consists in assigning 
arbitrary values to two independent variables from which an exact 
equilibrium-composition total pressure and quantity of oxygen are 
computed. By a process of successive approximations, the values of 
the two independent variables are so adjusted that the calculated 
total pressure and quantity of oxygen approach the desired values. 
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Charts are presented that aid in the choice of the initial estimates 
and rapidly effect their adjustment. Additional charts are included 
that may he used to calculate the partial pressures of the products 
of reaction. 

The compounds that are to react are represented by the equivalent 
formula C y E x OyS z where w, x, y, and z are proportional to the 

total number of carbon, hydrogen, oxygen, and nitrogen atoms con- 
tained in the initial compounds. For convenience, the symbols are 
listed in appendix A. . 

The reaction under consideration can then be expressed in the 
following form: 

A <W/ Z ) ^ P H 2 0 + p Hj + p 0H * % + + %(> + % + p C0 2 + %> 

+ Po 2 + Po (D 

where A is the number of formula weights of the equivalent for- 
mula C w H x O y N z and the volume of the products of reaction V is 

numerically equal to the product of the gas constant and the tempera- 
ture ET so that the total number of moles n of gaseous equilib- 
rium products is numerically equal to the total pressure P. The 
partial pressure p^ of each gaseous product is similarly equal to 

the number of moles n^ of the product i, for example, 

\o s % 0 * 

It is shown in appendix B that the partial pressure of each of 
the products can be computed when the value of the ratio of the 
partial pressure of water to the partial pressure of hydrogen, here- 
inafter designated the water-hydrogen ratio, ¥ = Ph^o/ph^ = KqPcOg/PCO, 

and the value of A are known. The equilibrium constant K^_ is 

defined by equation (19), appendix B. The partial-pressure calcula- 
tions- can be arithmetically made by the use of equations given in 
appendix B, or if the required accuracy need not exceed 2 percent, 
the values may be found from the calculation charts shown in fig- 
ure 1. The operations required for the use of figure 1 are dis- 
cussed in a subsequent section where the details of its use are 
illustrated by an example. 

In actual calculation, values of A and "P are initially 
estimated and a corresponding equilibrium composition is computed. 
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From the trial equilibrium composition, the corresponding total pres- 
sure P' and oxygen quantity Ay' are computed from the relations 

P' = PHgO + p Hg + p OE + PH + PNg + PNO + PN + P C0 2 + P C0 + POg + p 0 

( 2 ) 

and 

Ay ' ■ "HjO + P OH 2p C0 2 + p co + 2p 0 2 + p 0 <3) 

where primes denote specific values of the variable resulting from 
trial calculations. 

If the total pressure and the oxygen content resulting from a 
trial calculation are sufficiently close to the desired values, the 
problem is complete; otherwise the values of A and 9 must be 
adjusted. 


ESTIMATION AND ADJUSTMENT OF A AND 9 


Estimation. - Instead of directly estimating A, estimation 
of the dissociation factor AZ/P^ is more expedient where 


Z = w + x/2 + z/2 


and 


P 


1 


= P 



“ PO 


(4) 


(5) 


The value of AZ/P^_ is always between l/2 and 1 as shown in 
appendix B. 

The oxygen ratio r is defined by the relation 

_ 2w + x/2 ( 6 ) 

y 

The quantity r is used to indicate mixture ratio and is the ratio 
of oxygen required for stoichiometric combustion to that supplied, 
therefore r is equal to 1 for a stoichiometric mixture and 
approaches zero as the proportion of oxygen is increased. 
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Values of ¥ and AZ/P-^ are plotted in figure 2 as a function 

of r and temperature T for a pressure of 20.4 atmospheres for 
mixtures of oxygen with carbon, hydrogen, and hydrocarbon that has 
a hydrogen- carbon ratio x/w of 3. The effect of nitrogen on the 
values of ¥ and Al/P^ is sufficiently small that it may be neg- 
lected in making the first estimates of ¥ and AZ/P^. 

The desired values of r and Z are computed from the specified 
values of w, x, y, and z. Values of ¥ and AZ/P-^ are then 

estimated from the part of figure 2 that corresponds most closely to 
the x/w ratio of the reactant. (A considerable latitude in the 

choice of the initially estimated value of ¥ is permissible because 
the recommended process of adjustment rapidly corrects for initial 
errors. ) 


The partial pressures Pq and Pq are then computed by means 

2 

of equations (36) and (37) or by the use of figure 1(a). By use of 
equation (5), P 1 is then computed and A is computed from the 

expression 


A - (AZ/P-l) Pj/Z (7) 

Adjustment ♦ - The initially estimated values of ¥ and 
AZ/P-l will, in general, require adjustment. At the end of a trial 

calculation, a more exact value for each quantity can be graphically 
obtained with aid of figure 2 or a similar plot. 

The quantities AZ/P^’ and ¥ can be plotted on the appropriate 
part of figure 2 against the quantity r ' 

where 

r* = A (2 w + x/2) (8) 


V - p ’ - " *0 

and the quantities P* and Ay 1 are computed from equations (2) 
and (3). A curve, similar to the existing curves, is drawn through 
the plotted points. New values of ¥ and AZ/P^ corresponding 

to the desired value of r are read from the newly sketched curve. 
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The process of obtaining and adjusting values for the estimated 
quantities is illustrated in the section EXAMPLE OF EQUILIBRIUM - 
COMPOSITION CALCULATION. 

Properties of quantities ¥ and AZ/P^. - The values of ¥ and 
A2 /P x for the three groups of reactants , carbon-, hydrogen-, and 

hydrocarbon- oxygen (hydrocarbon has an x/v ratio of 3), at a tem- 
perature of 3500° K and a pressure of 20.4 atmospheres are shown in 
figure 2(d). 

For values of r less than 1.7, the value of ¥ is principally 
a function of r and to a first approximation is independent of ¥ and 
the x/v ratio. For values of r greater than 1.7, the value of ¥ 
and the shape of the ¥- isotherm critically depend on the x/w 
ratio of the reactants. 

The values of Al/P^ and the shapes of the isotherms of AE/P ^ 

are nearly independent of the x/v ratio unless the proportion of 
carbon becomes large. For a reaction involving only carbon and 
oxygen, the quantity AE/P-l always equals 1. The values of ¥ and 

AE/P^ for other x/v ratios may be estimated by visual interpolation 

between figures 2(a), 2(b), and 2(c) or, if desired, approximate curves 
of ¥ and AE/P! may be computed for other hydrogen-carbon ratios 
by the method that is discussed in appendix B. 

An asymptotic property of the quantity ¥ may be used to 
estimate ¥ and AZ/P^ for total pressures other than 20.4 atmos- 
pheres. As r is reduced (proportion of oxygen increased) the 
proportion of diatomic and monatomic oxygen in the products increases; 
at the limit of r * 0 only these species are present; hence 

P = Po 2 + Po (10 > 

For a given value of temperature and pressure, ¥ can be determined 
from the relation 

P = Kg'? 2 + K 7 ¥ (11) 

This value of ¥ can be used to define a pseudo temperature T^ 
for which the pressure equals 20.4 atmospheres, that is 

20.4 . K 6jI / + K 7>p ? 


( 12 ) 
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where K g and K 7 ^ are the values of K g and K^, respectively, 

at the required pseudo temperature. Experience has shown that the 
pseudo temperature computed at r = 0 may he used for other values 
of r to obtain nearly correct initial estimates of both ¥ and 
AZ/P]_ at pressures other than 20.4 atmospheres. 

A plot of Tp against T for a range of pressures is shown 

in figure 3. In order to obtain estimates of ¥ and AZ/P-^ at 

pressures other than 20.4 atmospheres, the value of corresponding 

to the desired assigned temperature and pressure is found from fig- 
ure 3. This Tp is used in figure 2 instead of T to obtain the 

first estimates of ¥ and AZ/P-^. 

EXAMPLE OF EQUILIBRIUM-COMPOSITION CALCULATION 

The procedure used to compute equilibrium composition is illus- 
trated by the calculation of the composition of the products of 
reaction of ethylenediamine ( C g ( NHg ) g ) and oxygen (O2) at the stoi- 
chiometric ratio, a temperature of 4000° K, and a pressure of 10 atmos- 
pheres. The reaction is 


C 2 H 4 (NHg) 2 + 4 0g — ^-products 


By counting the number of atoms of each element in the equa- 
tion, the following values of w, x, y, and z are obtained: 

w * 2, x = 8, y = 8, z=2 

From these values and equation (6), the desired r Is 1 (that is, 
the stoichiometric ratio): 


r = 2* + x/2 _ 

y 

and from equation (4), Z is 7: 

S = w + x/2 + z/2 = 7 
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For the assigned pressure of 10 atmospheres and the assigned 
temperature of 4000° K, a T p of 4200° K is found from figure 3. 

For the desired r of 1 and a T p of 4200° K, the initial estimates 
of ¥ and AZ/P-l are found from figure 2(c) to he 1.6 and 0.8, 
respectively. The estimated value of AS/P^ is very close to the 

correct value; however, to illustrate the method of adjustment and 
the rapidity of its convergence, a value of A^/Pp = 0.84 was 

deliberately chosen for the computation. In order to illustrate 
this example, an enlarged section of figure 2(c) has been reproduced 
as figure 4 and the points required for this problem are plotted. 

From equations (36) and (37) of appendix B and the values of 
corresponding to the estimated value of ¥ of 1.6 are found to be 

Pq = 0.8082 

2 

Pq = 1.3938 


and IL, in table I for 4000° K, the values of p Q and p^ 


and by equations (5) and (7) 


P, = P - p_ - p n = 10 - 0.8082 - 1.3938 = 7.7980 
1 Uo u 


and 

A = (AZ/P 1 )(P 1 A) * 1^84)^7/798} „ 0 .9358 

By the use of equations (39), (40), and (42) to (50) of appendix B 
and the equilibrium constants of table I corresponding to 4000° K 
the remaining partial pressures are found to he: 


p =3 

= 0.8782 


= 0.7907 

f C0 2 

= 0.2893 

P 0H 

= 1.4236 

P N0 

= 0.2409 

p co 

= 1.5823 

PXT 

= 1.4965 

P» 

= 0.0493 




P V = 


1.4051 
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From these partial pressures and equations (8), (2), and (9), 
the values of r*, P', and AZ/Pj • corresponding to the estimated 

values of A and ¥ are found to be 

r t e A(2w + x/2) 

Ay* 

= 0.9085 


where 


Ay ' “ ^HgO + P 0H + P N0 + 2p C0 2 + P C0 + 2p 0 2 + p 0 
= 8.2407 


and 

P'=p„ +p +p +P+P +P +P+P +P +P +P 
HgO Hg OH H N 2 NO N C0 2 CO 0 2 0 

= 10.3579 atmospheres 

* 10.3579 - 0.8082 - 1.3938 
= 8.1559 


therefore 


AZ/P, , . (0.9358) (7) 

1 8.1559 

= 0.8032 

If the values of r' and P' do not agree closely enough with 
the desired values, the quantities ¥ and A must be reestimated. 
The values of ¥ and AZ/P^ ' can be plotted in figure 2(c) against 

r’ and are again shown in figure 4. Segments of curves similar 
in shape to the existing curves are sketched through the plotted 
points; new values for ? of 1.4 and for AZ/P-^ of 0.81 are read 

from the newly sketched curves at r = 1 (the desired value). 
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The second trial calculation is based upon the newly estimated 


values 

of ¥ and 

AZ/P-l with 

the result that 



P °2 

= 0.6188 

p„ = 0.9728 
“2 

P »2 

= 0.8127 

P C0 2 

= 0.2605 

p o 

= 1.2195 

P 0H = 1.3111 

P N0 

= 0.2137 

O 

O 

= 1.6283 

P 1 

= 8.1617 

P H = 1.5751 

% 

= 0.0500 



A 

= 0.9444 







whereby 
r' = 1.0083 

P' = 10.0244 atmospheres 
AZ/P-l* = 0.8076 


These values are in close agreement with the desired values of r 
and P. 

The detailed calculations and the results of this example are 
shown in table II. 


GRAPHICAL CALCULATION OF PARTIAL PRESSURES 

The numerical calculation of partial pressures involves tedious 
and time-consuming operations that include reference to tables of 
equilibrium constants. As an aid to the computation, the solutions 
of the equilibrium equations have been plotted in figure 1. 

The procedure required to obtain the partial pressure of each 
of the products of reaction from the charts is illustrated for the 
preceding example by the guide lines shown on each part of figure 1. 
The conditions for the example are 

C ? H 4 (NH2) 2 + 4 0 2 — > products 

w = 2, x = 8, y = 8, z = 2 
r s 1 ,Sb 7, P ■ 10 atmospheres, T = 4000° K 


MCA TN No. 1653 


11 


As in the preceding example, the first estimates of ¥ and A2/P^ 
are 1.6 and 0.84, respectively. 

The partial pressures Pg and pg are plotted against ¥ 

2 

in figure 1(a) for various temperatures. For ¥ * 1.6 and T * 4000° K 
p Q m 0.81 and Pg - 1.39. 

2 

Whereby P x ■ 10 - 0.81 - 1.39 - 7.80 

A « 0.84 7 -*-§QQ - 0.936 
7 

Jx _ (0.936X8). = 2 8 

¥+1 1 . 6+1 


Az « (0.936)(2) = 1.872 


Aw « ( 0 . 936 ) ( 2 ) = 1.872 


A plot of ¥ against a, which is defined by equation (39), 
for various temperatures and a plot of Pg against a for various 

Lt 

values of Ax/¥+l are given in figure 1(b). For ¥ * 1.6 and 
T = 4000° K, a = 1.20; for a = 1.20 and Ax/¥+l - 2.88, 

p_ «= 0.88. 

*2 


Figure 1(c) is composed of two parts. The constant tempera- 
ture lines in the upper left section are plots of pjj against 
p H ^ for various temperatures. For p^ « 0.88 and T *= 4000° K 

(read from the scale marked temperature T for atomic hydrogen H) , 
p H » 1.50. By the use of a set of reference lines that correspond to 


constant values of 



and the second set of constant tempera- 


ture lines (in the lower right section) p QE is computed. For 


¥ m 1.6 and pg^ = 0.88, the reference line q is found. At the 

intersection of q and the temperature line T = 4000° K (read from 
the scale marked temperature T for hydroxyl radical OH), Pgg * 1.42 


A plot of ¥ against b, which is defined by equation (42), for 
various temperatures and a plot of p against b for various 


values of 
T 


4000° K 


»2 

Az are shown in figure 1(d). 
b = 0.326. For b = 0.326 


When ¥ = 1.6 and 
and Az • 1.872, p^ 


0.79. 
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The constant temperature lines In the upper left section of 
figure 1(e), which Is composed of two parts, are plots of pjj against 

Pjj for various temperatures. For p^ = 0.79 and T * 4000° K 
2 2 

(read from the scale marked temperature T for atomic nitrogen N), 
p^ = 0.049. By the use of a set of reference lines that correspond 


to constant values of ¥ 




and a second set of constant tempera- 


ture lines (in the lower right section) p^Q is computed. For 

¥ = 1.6 and p^. = 0.79, the reference line q is found. At the 
w 2 

intersection of q and the temperature line T = 4000° K ( read 
from the scale marked temperature T for nitric oxide NO), 

Pno = °* 241, 

A plot of p^Q /Pqq against ¥ for various temperatures and 
2 

plots of p co ^ against P nn _/Pp n for various values of Aw when 


TO*' *C0 


Pco /p co _ 1 and p co plotted against Pqq /vqq for various values 
2 2 

of Aw when p AA /p ^ 1 is presented in figure 1(f). This chart 
OO 2 C0 

is so constructed that only the smaller of the two quantities 

p or p AA is determined. The larger of either of these quanti- 

COg co 

ties is obtained by subtraction from Aw. For ? « 1.6 and T * 4000° K, 


/p = 0.183 and for Aw = 1.872, p 


•C0r/ ^C0 


CO. 


0.290 and 


P co « 1.872 - 0.290 


1.582. 


The quantity p_ is equal to f times p„ ; therefore 
HgO h 2 

= (1.6) (0.88) = 1.408. From these values of partial pressures, 
Ay' = 8.241, r' « 0.9086, and P-^ « 10.360. 



New estimation of 5* and A are obtained as previously 
explained and a second trial calculation is made if the values of 
P' and r* are to be brought closer to the desired values of 10 and 
1, respectively. The results of two trial calculations made by the 
graphical method have been tabulated in table II. 

Numerical checks on calculations. - A check on the accuracy of 
the chart calculations can be made by substituting the partial pres- 
sures obtained from a trial calculation in equations (32) and (34) 
to determine Ax and An. The values obtained for this example are 


NACA TN No. 1653 


13 


Ax * (2) (1.408) + (2) (0.88) + 1.42 + 1.50 = 7.496 

Az = (2)(0.79) + 0.241 + 0.049 « 1.870 

These values axe to "be compared with Ax = (0.936)(8) « 7.488 and 
Az * 1.872. 

When the values Ax and Az do not check exactly, improvement 
of the accuracy of the calculations is possible by use of equations 
for the mass balance. 


If, as will usually be the case, the slam of p^ and p^ is 

less than -^Ax and the sum of p and p„ is less than -iAz, 

2 NO N 2 ' 

the following equations can be used to compute better values for 
p , p , and p 

V %* %0 



Az ~ %> ~ % 
2 


(13) 


1.872 - 0.241 - 0.049 
2 


= 0.791 


when 


\o> % 


p * ~ ^ - p 0H - p H 

% 2 (? + 1 ) 


(14) 


7.488 - 1.42 - 1.50 

( 2 ) ( 1 . 6 ) + 1 



= 0.878 

Ax - p 0H - p H ^ 

2 Pfl 2 


( 15 ) 


s 7_i^§§_I-JLt£§_I_ii^2 — 0 878 
2 


1.406 
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and when 





Ax - Pqh - Ph ? 

2(f + 1) 


(16) 



** " P 0H " P H _ 

i Pa 2° 


(17) 


Flight Propulsion Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, April 1, 1948. 
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APPENDIX A 
SYMBOLS 

A 

number of formula weights of the equivalent for- 
muia 

AZ/P 1 

disBOciation factor, moles/atm 

a 

K 9 * * K 10 

— ^ ^ term of quadratic equation 

b 

K Q + WKj term of quadratic equation 

WA 

equivalent formula for compounds that are to 
react 

... 

equilibrium constants 

n 

total number of moles of products of reaction 

n H 2 0' n ^2 > * * 

. number of moles of water, hydrogen, . . . , 
respectively 

P 

total pressure of products of reaction 

P 1 

p - p 0 2 - p o 

V' V ' * 

. partial pressure of water, hydrogen, . . . , 
respectively 

R 

gas constant, 1.987, cal/ (mole) (°K) 

r 

. . 2w + x/2 

oxygen ratio, 

3 

T 

temperature, °K 

T 

P 

pseudo temperature, °K 

V 

volume of products of reaction, cal/atm 

v, w, x, y, z 

total number of atoms of lithium, carbon, hydrogen, 
oxygen, and nitrogen, respectively, in initial 
compounds. 
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Z v + x/2 + z/2 

f water- hydrogen ratio. 

Primes are used to indicate specific values of the variable 
resulting from trial calculations. 



G 
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APPENDIX B 


DERIVATION OF EQUATIONS 


In this section it is shown how the equations for the conservation 
of mass and the equilibrium equations are combined to form the work- 
ing equations required for the computation of an equilibrium com- 
position. The range for the quantity A2/P 1 and expressions used 

for plotting approximate and exact hydrocarbon curves, similar to 
those of figure 2, are shown. 


The following dissociation reactions and the corresponding equilibrium 
equations were considered for the derivation. These equations 
together with the values of the equilibrium constants were taken 
from reference 3. 


Basic Equations 


By definition 



(18) 


C0 2 + Hg^^CO + HgO 



(19) 



i^=±NO + Hg 



%0 


( 20 ) 



*6 * f 2 
HgO^O + Hg 


( 21 ) 



( 22 ) 
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•=N f > N 
2 2* 


K, 


TI 


8 


1 = 2 ^ 


K 9- 


r H 




(23) 


(24) 


*2^*2 + 


OH 


*10 ” 


•OH 






(25) 


The lav of conservation of mass yields the following equations: 
A ( C v^ic 0 y W z^ — * d f o + ^2 + “OH + n H + + "NO + °N + “COg + “CO 


+ no 2 + no 

Ax = 2n s ^ 0 + 2ng 2 + Dqh + n H 
Ay " n HgO + "OH + "NO + 2n C0 2 + "CO + 2n 0 2 + "O 
Az = 2 "n 2 + "NO + "N 
Aw * "COg + "CO 


(26) 

(27) 

(28) 

(29) 

(30) 


Equations (18) to (25) are in terms of partial pressures whereas 
equations (26) to (30) are in terms of moles. When a volume is so 
chosen for the reaction that 


m « i 


and the perfect-gas law is assumed, that is 

PV = nRT 

P BT . 
n * V 1 

and 


( 31 ) 


n 
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Similarly 


P H 0 = "h.O 
2 2 


therefore equations (27) to (30) may he written 

V + 2PH 2 


Ax = 2p LQ + 2 Ph _ t, P QH + P H 


Ay = P H 2 0 + P 0H + P N0 + 2P C0 2 + P C0 + 2P 0 2 + P 0 


AZ " 2p H 2 + P N0 + P N 


Aw = P C0 2 + P C0 


(32) 

(33) 

(34) 

(35) 


Working Equations 

The equations shown in the preceding section are combined to 
form the relations used for calculating the partial pressures. The 
quantities p. and p. are obtained from equations (2l) and (22). 


-*% 

(36) 

tt 

M 

(37) 


Equations (18), (24), (25), and (32) are combined to form the 
following quadratic equation for P g^ in terms of ?, Kg, K^q, 

and A: 


K 9 +fK 10 . 

2Ptt + i/P 






Ax 


Ho 


¥ + 1 


= 0 


? + 1 


(30) 
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if 

a = K 9 + ^0 (39) 

? + 1 




(40) 


Equations (20), (23), and (34) are combined to form the follow- 
ing quadratic equation for Pjj in terms of W, K 3 , Kg, and A 


2 % 2 + ^8 + ' PK 3 ) f% 2 - Al '° (41) 


if 



(42) 

(43) 


For given values of x, t, and temperature, and by estimating 
values for A and ¥, p^ and p^ can be computed. 

The remaining partial pressures can then be computed. Equa- 
tions (18), (20), and (23) to (25) are rewritten: 


v " ?ph 2 

(44) 

P 0H = fE 10 ^ P Hg 

(45) 

P H ' K 9 f % 

(46) 

%0 ‘ tK 3 lJ\ 

(47) 




( 48 ) 
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By combining equations (19) and. (35) 


„ _ Av 

C0 3 ‘ 775 


p co * A “ - p co. 


or 


Aw 


CO . ¥ 

1+ 7 


p co 2 ' Av - p co 


Range of Values of AZ/P-^ 

From equations (2) and (5) 

P=p +p +p + p + p +p +p+p +p 

1 \o \ y 0E P H ^N 2 P N0 P N P C0 2 P C0 

and from equations (4), (32), (34), and (35) 

AS = Aw + Ax/2 + Az/2 

x r x-n x-n ^ P 0H + PH JT , ^PNO + PN 

P C0 2 + P C0 + P Hg + % 2 + “ 

By division of equation (52) by (51) and by simplification 

P 0H + P H + ^0 + % 


A^ = 1 - 


ZP, 


The maximum value of AZ/P occurs where p + p + p + p 

1 OH H NO 

and the minimum value occurs when p QH + Pg + Pjj 0 + Pn “ F l J 


(49) 

(50) 

(49a) 

(50a) 

(51) 

(52) 

(53) 
hence 
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1 ^ az/p-l = | 

AE/P! Is a dissociation factor, the value of which is less than 1 

when molecules containing one atom of hydrogen or nitrogen are 
present. This factor is used to estimate A. 


Approximate Interpolation Formulas for 
Hydrocarbon-Oxygen Reactions 

Approximate curves of ? and AE/P^ as a function of r can 
be constructed for any hydrocarbon-oxygen reaction by the following 
method: 

For a given value of ¥ and temperature, r^ and rjj are 
read from figures 2(a) and 2(b), respectively, and (AE/P-jJg is 
read from figure 2(b) at r^. Values of rgg and (AE/P^Jqj for 
the hydrocarbon are computed from the approximate equations. 



The quantities ¥ and (AZ/Pj^qjj are plotted against Tqjj to 
form a family of curves similar to those of figure 2(c). 

Equations (54) and (55) are exact for the limiting case of low 
temperature and high pressure when p QH and p 0 are equal to zero; 

equation (54) is exact for the limiting case of high temperature and 
low pressure when PHgO 8X1(1 Pflg are to zero. 

Except for special cases, small errors exist in both equa- 
tions (54) and (55). The order of magnitude of these errors is 
indicated in the following table for T = 4000° K, and P *= 20.4 atmos- 
pheres. 
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x/w 

¥ 

lea 


Exact 

value 

Approximate 

value 

Error 

(percent) 

1 

0.1 

2.3425 

2.3454 

0.12 

1 

.3 

2.1103 

2.1183 

.38 

1 

.4 

.6499 

.6718 

3.37 

4 

.1 

3.2648 

3.2893 

.14 



(az/p 1 ) ch 


1 

0.1 

0.91584 

0.91352 

0.25 

1 

.3 

. 91396 

.91152 

.27 

1 

.4 

.86366 

.85657 

.82 

4 

.1 

.87814 

.87770 

.05 


Exact Equations for Plotting Hydro cart on Curves. 

The following exact equations may te used to construct plots 
of A2 /P j and ¥ against r similar to those in figure 2(c) for 

the carb on- hydrogen- oxygen system: 


r CH 


2¥ + kA 
+ k i. 


2w + x/2 

x? 

+ i) (2 + “) 



+ K io 


+ 


1*% ) 



(56) 

(57) 


where 
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Equation (58) may be used to simplify the calculation of 
equilibrium compositions in the carbon- hydrogen- oxygen system. 

When equation (58) is used, the total pressure is always exact and 
it is therefore unnecessary to estimate AZ/P^j only f need be 
estimated. 
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APPENDIX C 

CALCULATION OF EQUILIBRIUM COMPOSITION WITH ELEMENTS IN 
ADDITION TO CARBON, HYDROGEN, OXYGEN, AND NITROGEN 


The method presented for the solution of the carbon-hydrogen- 
oxygen-nitrogen system often can be used to aid in the solution of 
problems involving other elements by dividing the reaction into two 
parts. For example, if lithium Li is added to the carbon-hydrogen- 
oxygen-nitrogen mixture, the reaction can be expressed by the equation 

Li v C w H x O y N z -^ products 
where v is the number of atoms of lithium. 


Let the first part of the reaction be 

C w H x O y X N Z — ^-products of C, H, N, and 0 

and the second part will be 


LLyOy 2 — > products of Li and 0 

In order to make this division of the problem, the following require- 
ments must be met: 


1. The pressure of 



and 


Pq must be the same in both 


parts of the reaction and must be considered only once in the com- 
putation of the total quantity of oxygen or the total pressure. 


2. The formation of compounds other than oxides must be con- 
sidered and, if present, they must be correctly accounted for 
(for example, LiH). 

The products of the first part of the reaction can be found 
by means of figure 1 and the products of the second part are determined 
by some other method. 

If it may be assumed that no reaction except the formation 
of the oxide is present, the method becomes fairly easy. An example 
is subsequently presented to illustrate the procedure for calcula- 
tion. For this prupose, a dissociation reaction has arbitrarily 
been assigned: 
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LigO -f < — — HgO + 2 Li 


the dissociation constant has been assigned an arbitrary value 
because no data are available, therefore 


ty assumption for a temperature of 4400° K. All compounds and ele- 
ments are assumed to be gases. These assumptions probably are not 
valid for the lithium reaction and should not "be used for any purpose 
other than illustration. 

Let the desired reaction he: 


v = 2, i = 4, y = 3, T = 4400° K, P = 20.4 atmospheres 
Because P = n 



K = — = 0.1 


(59) 


P H 2 ^2° 



for 


2 PLi 2 0 + PLi = Av 


(60) 


When equations (18), (59), and (60) are combined 



(61) 


let 



(62) 


and 



( 63 ) 
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Av - 


^Li 


PLigO *= “ _ 

For this reaction r would, be expressed 

v + x 


r = 


2y 

2 + 4 


(64) 


= 1 


and E would be 

E = v/2 + x/2 = 2/2 + 4/2-3 

where v/2 appears because two atoms of Li are required to form 
the normal oxide LigO. In the absence of a plot of ¥ and 

AE/Pj_ for lithium, hydrogen, and oxygen, figure 2(b) may be used 
from which 


¥ =0.87 and AE/P 1 = 0.715 
and from figure 1(a) 


p = 0.96 
°2 

p Q = 3.05 

whence 

P = 20.4 - 0.96 - 3.05 = 16.39 

A - (AE/P X ) Pj/E = 0.715 i 6 .-, 5 .3= 3.91 

3 

Ax/(¥ + 1) = <yiK4 i = 8.36 


From figure 1 
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2.25 

(2.25) (0.87) 


1.96 


P 0H = 2.72 
P H = 4.51 


From equation (62) 


0.1 

(4) (o.87 ) 


0.0287 


and 


AvK (5.91)(2)(0.1) 
2¥ " (2) (0.87) 


0.449 


and from equation (63) 

= -0.0287 + V(0*0287) 2 + 0.449 = 0.642 

From equation (64) 

_ (3.91) (2) - 0.642 = 3>59 

Li 0 2 

2 

P 1 = 2.25 + 1.96 + 2.72 + 4.51 + 0.64 + 3.59 = 15.67 

P = ? 1 + 0.96 + 3.05 = 19.68 

Ay = (2)(0.96) + 3.05 + 1.96 + 2.72 + 3.59 = 13.24 


A(v + x) (3.9l)(6) 
r = 2Ay " (2) (13.24) 


0.886 


AS/?! 


(3,91) (3) 
15.67 


0.749 


The composition thus obtained is correct for a pressure of 19.68 atmos- 
pheres and for an oxygen ratio of 0.886. A closer approximation to 
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the desired values of r - 1 and P . 20.4 atmospheres can be 
obtained by the method of adjustment of A and l deeorlbed for 
the calculation of composition for the carbon- hydrogen-oxygen- 
nitrogen system. 
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TABLE I - EQUILIBRIUM CONSTANTS 
[Data obtained from reference 3. ¥ ■ PH 20 /PH 2 ] 


Temper- 

ature 

(°K) 


P C0 y 

Pn° 

P °2 

K 6-*f 

K-^ 


lJI. 

P ° H 

2 

Vpn^ 

% m; 

9 mi 

10 *'P Vp£ 

4.782 

5.736xl0 -6 

8.139x10" 8 

2.006xl0 -7 

8. 874x10 ~ 7 

1.623xl0 -3 

1.636x10“ 

5*149 

1.531x10*5 

3.425x10* 7 

8.847xl0 -7 

2.528x10”® 

3.109x10”® 

3.682x10“ 

5*420 

3.785x10”® 

1.307x10- 6 

3.345x10”® 

6.544x10”® 

5.620x10”® 

8.031x10“ 

5.853 

8.373xl0“ 5 

4.173xl0’ 6 

1.093x10”® 

1.562x10”® 

9.662x10”® 

1.561x10“ 

6*140 

1.771x10"^ 

1.260x10*5 

3.307x10”® 

3.469x10”® 

1.588x10-2 

2.939x10“ 

6.440 

3.512xl0" 4 

3.452x10*5 

9.121x10”® 

7.231x10”® 

2.514xl0”2 

5.230x10“ 

6 668 

6.644xl0* 4 

8.848x10*5 

2.337xl0” 4 

1.426xl0” 4 

3.835x10-2 

8.954x10* 

6*954 

1.188x10*5 

2.080x10*4 

5.546x10”* 

2.672xl0” 4 

5.674x10-2 

1.461x10* 

7.185 

2.047xl0* 3 

4.639x10*4 

1.243x10”® 

4.795X10” 4 

8.173x10”2 

2.308x10* 

7.376 

3.414x10*5 

9.883x10*4 

2.646x10”® 

8.267xl0” 4 

.1148 

3.550x10' 

7.582 

5.479x10*5 

1.983x10*5 

5.336x10”® 

1.373x10”® 

.1579 

5.290x10 

7.748 

8.525x10"^ 

3.803x10*5 

1.027x10“2 

2.210x10”® 

.2125 

7.677x10 

7.929 

1 .294x10*2 

7.055x10*5 

1.904x10”2 

3.452x10”® 

. 2813 

.1087 

8.078 

1.859x10“;? 

1.248x10*2 

3.389x10”2 

5.253x10”® 

.3654 

.1513 

8.192 

2.758x10’;? 

2.157x10*2 

5.857xl0”2 

7.798x10”® 

.4682 

.2060 

8.304 

3.894xl0“ 2 

3.583x10*21 

9.778xl0”2 

1.132x10-2 

.5910 

.2763 

8.449 

5.404xl0’ 2 

5.821x10*2 

.1591 

1.611xl0 -2 

.7367 

.3614 

8.535 

7.355x10*2 

9.166x10*2 

.2516 

2.248x10”2 

.9084 

.4691 

8.664 

9.821x10*2 

.1412 

.3876 

3.086xl0”2 

1.107 

.5984 

8.712 

.1300 

.2134 

.5878 

4.166xl0”2 

1.336 

.7589 

8.752 

.1693 

.3157 

.8711 

5.544x10”2 

1.597 

.9495 

8.819 

.2162 

.4526 

1.258 

7.258x10”2 

1.893 

1.164 

8.884 

.2760 

.6519 

1.806 

9.423x10-2 

2.226 

1.428 

8.898 

.3463 

.9139 

2.537 

.1206 

2.598 

1.730 

8.929 

.4302 

1.263 

3.508 

.1528 

3.012 

2.089 

8.952 

.5291 

1.718 

4.778 

.1916 

3.470 

2.478 

8.974 

.6451 

2.307 

6.431 

.2379 

3.969 

2.923 

8.984 

.7799 

3.068 

8.545 

.2927 

4.520 

3 . 432 

8.994 

.9344 

4.015 

11.21 

.3570 

5.118 

3.992 

9.003 

1.115 

5.228 

14.59 

.4321 

5.764 

4.637 

9.012 

1.317 

6.693 

18.77 

.5193 

6.461 

5.317 


2000 

2100 

2200 

2300 

2400 

2500 

2600 

2700 

2800 

2900 

3000 

3100 

3200 

3300 

3400 

3500 

3600 

3700 

3800 

3900 

4000 

4100 

4200 

4300 

4400 

4500 

4600 

4700 

4800 

4900 

5000 


.-4 


-3 

-2 


KM 

O 


NACA TN No. 1653 


NACA TN No. 1653 


31 


TABLE II - RESULTS OF SAMPLE CALCULATION OF THE REACTION 
c 2 h 4 (nh 2 ) 2 + 40 2 — PRODUCTS 

2, X = 8, y = 8, z = 2, Z = 7, r = l.OO, T = 4000° K, P = 10, 
8.752, K 3 = 0.1693, Kg = 0.3157, K? = 0.8711, Kg = 0.0554, 


K 9 = 1.507, K 


10 


0. 9495] 


Quan tl ty 

Ari thmetical 


Graphical 


Source 

First 

trial 

Second 

trial 

Source 

First 

trial 

Second 

trial 

9 

Figure 2(c) 

1.6 

1.4 

Figure 2(c) 

1.6 

1.4 

A2/P! 

Figure 2(c) 

.84 

.81 

Figure 2(c) 

.84 

.81 

p o 2 

Equation (36) 

.8082 

.6188 

Figure 1(a) 

.81 

.62 

P 0 

Equation (37) 

1.3938 

1.2195 

Figure 1(a) 

1.39 

1.22 

pi 

Equation (5) 

7.7980 

8.1617 

Equation (5) 

7.80 

8.16 

A 

Equa ti on ( 7 ) 

.9358 

.9444 

Equation (7) 

.936 

.944 

Ax/( •F+l) 
a 

(a) 

Equation (39) 

2.8793 

1.1985 

3.1480 

1.2193 

(a) 

2.88 

3. 146 

Ph 2 

Equation (40) 

.8782 

.9728 

Figure 1(b) 

.88 

.97 

P 0H 

Equation (45) 

1.4236 

1.3111 

Figure 1(c) 

1.42 

1.31 

% 

Equation (46) 

1.4965 

1.5751 

FI gure 1(c) 

1.50 

1.57 

Ph 2 o 

Equation (44) 

1.4051 

1.3619 

Equation (44) 

1.408 

1.358 

Az 

b 

(a) 

Equation (42) 

1.8716 

.3263 

1.8888 

.2925 

(a) 

1.872 

1.888 

% 

Equation (43) 

.7907 

.8127 

Figure 1(d) 

.79 

.81 

PnO 

Equation (47) 

.2409 

.2137 

Figure 1(e) 

.241 

.214 

P N 

Equation (48) 

.0493 

.0500 

Figure 1(e) 

.049 

.05 

Aw 

(a) 

1.8716 

1.8888 

(a) 

1.872 

1.888 

P C0 2 

Equation (49) 

.2893 

.2605 

Figure 1(f) 

.290 

.260 

p co 

Equation (50) 

1.5823 

1.6254 

Equation (50) 

1.582 

1.628 

P» 

Equation (2) 

10.3579 

10.0244 

Equation (2) 

10.360 

10.01 

Ay' 

Equation (3) 

8.2407 

7.4931 

Equation (3) 

8.241 

7.49 

r' 

P 1 , 
ax/p; 

Equation (8) 
Equation (9) 

ax/p; 

.9085 

8.1559 

.8032 

1.0083 

Equation (8) 
Equation (9) 

AX/P{ 

.9086 
8. 160 

.803 

1.008 


a 


Numerical calculation. 


Partial pressure of oxygen, p 02> and atomic oxygen, p 0 , 
. atmospheres 
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(a) Partial pressures of oxygen Pq 2 and atomic oxygen p Q . 

Figure 1. - Partial- pressure computation charts. ¥ « 17 " by 

20-inch print of this figure is attached.) 


G 
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atmospheres 

(b) Partial pressure of hydrogen pjj . Number of atoms of 

hydrogen Ax. 

Figure 1. - Continued. Partial-pressure computation charts. 

¥ - p h 2 o/ph 2 * 

(A 17- by 20-inch print of this figured attached.) 


of equation (39) 


Partial pressure of hydrogen, p H » atmospheres 
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Partial pressure of atomic hydrogen, Pjj, atmospheres 



(c) Partial pressures of hydroxyl radical p QH and atomic 

hydrogen p H » 

Figure 1. - Continued. Partial-pressure computation charts. 

¥ - Ph 2 o/ph 2 * 

(A 17- by 20-inch print of this figure is attached.) 


Partial pressure of hydroxyl radical, p QH , atmospheres 
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atmospheres 

(d) Partial pressure of nitrogen p™ . Number of atoms of 

nitrogen Az. 

Figure 1. - Continued. Partial- pressure computation charts. ¥ ■ Ph 2 o/Ph 2 * 
(A 17- by 20 -inch print of this figure is attached.) 


b of equation (42) 
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Partial pressure atomic nitrogen, p^ # atmospheres 



(c) Partial pressures of nitric oxide and atomic nitrogen 

p n . 

Figure 1. - Continued. Partial-pressure computation charts. ¥ ■ ^Hoo/PHo* 
(A 17- by 20-lnch print of this figure is attached.) 


Partial pressure of nitric oxide, p ^ q, atmospheres 


Partial pressure ratio of carbon dioxide and carbon 
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Partial pressure of carbon dioxide, p C q 2 , 81X1(1 carbon monoxide, 

P C q, atmosphere^ 



(f) Partial pressures of carbon dioxide Pc02 81X1(1 cart)on 

monoxide Pqq. Number of atoms of carbon Aw. 

Figure 1. - Concluded. Partial-pressure computation charts. ¥ ■ c/ p H * 

2 2 

(A 17- by 20 -inch print of this figure is attached.) 


Water-hydrogen ratio, $ Dissociation factor 

AE/P^, moles/atm 
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Temperature 

(°K) 






2000 to 50< 







JU 
































Figure 2. - Variation of water -hydrogen ratio ^ and dissociation 
factor AE/P^ with oxygen ratio r at a pressure of 20*4 

atmospheres, ¥ « PH20/PH2* (An by 17-inch print of this 

figure is attached.) 
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of 20.4 atmospheres. ¥ m Ph 2 o/PH 2 # (An 17-inch print 

of this figure is attached.) 


Water -hydrogen ratio, ^ Dissociation factor 

Ae/P^, moles/atm 
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Temperature 



.2 .3 .4 .6 .8 1 2 3 4 6 8 10 20 


Oxygen ratio, r 

(c) Hydrocarbon and oxygen with a hydro gen -to -car bon 
atom ratio of 3. 

figure 2. - Continued. Variation of water-hydrogen ratio ¥ and 
dissociation factor AE/P^ with oxygen ratio r at a pressure 

of 20.4 atmospheres. $ » ^Hoo/PHo* (An by 17-inch print 

c it 

of this figure is attached.) 
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.2 .3 .4 .6 .8 1 2 3 4 6 8 10 20 

Oxygen ratio, r 


(d) Carbon-, hydrogen-, and hydrocarbon-oxygen at 
temperature of 3500° K» 

Figure 2. - Concluded. Variation of water-hydrogen ratio ¥ and 
dissociation factor AE/P^ with oxygen ratio r at a pressure 

of 20.4 atmospheres. * Pho0^ p Ho* 

C* 
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Water-hydrogen ratio, i* Dissociation factor 
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Temperature 

(°K) 



Figure 4. - Application of figure 2(c) to example, ethylenediamine 
C 2 H 4 (NH 2 ) 2 reacting with oxygen 0 2 . ^ = Ph 2 o/Ph 2 * Pressure, 

10 atmospheres; temperature, 4000° K. 


NACA - Langley Field, Va. 




